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INTRODUCTION 


ARPSIM is a computer program developed to provide estimates of the terminal 
effectiveness of an Anti~-Radiation Projectile (ARP) fired against an air defense, 
vadar~enitting target. 


The primary objective of ARPSIM is to provide the user with a tool to para- 
metrically ascertain the sensitivity of the ARP to warhead, guidance, and fazing 
design changes. 


The ARPSIM model simulates single round terminal conditions frox the time 
when the ARP is flying a straight line trajectory at some fixed attack elevation 
in the vicinity of the target. Trajectories are determined from guidance errors 
distributed about a specified homing point. No further trajectory alterations 
are made. Fuzing points on the target are specified, and when fuzing conditions 
are satisfied, a burst poiat is established along the selected trajectory. The 
proximity of the burst poirt to the target determines the magnitude of kill pro~ 
babilities for blast, direci hit, and fragmentation effects. Separate blast 
kills for both the target body and radat antenna can be estimated. Fragmentation 
effects are based upon terminal effectiveness estimates generated by the full 
spray waterial lethal area (MAE) computer code (refs 1 and 2). 


The ARPSIM program is coded in FORTRAN for interactive use on the CDC 
6500/6600 in the INTERCOM mode. The user is prompted for data entry. Also, at 
the option of the user, an inpuc guide can be generated prior to each use. Frag~ 


mentation effects are estimated from data previously generated by the MAE program 
relative to conditional kill probabilities. Optionally, a function, Pi.(r), can 
be provided to specify fragmentation kill probability as a function of range. 
Comments are added throughout the FORTRAN code for better understanding and for 
development of future options for the code. 


A user guide, an example of a computer run, and a FORTRAN code listing are 
presented as appendixes A, B, and C. 


PROGRAM FLOW 


For each Monte Carlo sample, a simulation of the terminal characteristics of 
the ARP is made beginning at a time prior to fuzing during the ARP flight after 
final corrections to the trajectory have been made and when the remaining trajec- 
tory is linear at a fixed attack angle. The sequence of events for each simula- 
tion is: 


1. An attack angle is chosen which provides a straight line flight path 
with respect to a specified homing point. 


2. A trajectory is chosen based upon the guidance errors with respect to 
the homing point. 


3. A fuzing point along the chosen flight path is determined. 
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4. A burst point is established based on the type of fuze, direct hit 
potential and possible backup fuzing or ground impact prior to nominal fuzing. 


The proximity of the burst point to the target yields estimates of ki..l proba~ 
bility for direct hit, target body blast, radar blast, and fragmentation 
effects. The overall kill probability for each simulation is determined from the 
individual kill mechanism effects. This process is repeated for each simulation 
to provide the desired estimates of ARP terminal effectiveness. The above~de- 
scribed program flow is illustrated in figure 1. 


The following subsections briefly describe portions of the model in the 
approximate order in which they follow the program flow. 


Terminal Effects 


Terminal effects are measured in terms of direct hit, blast, and fragmerta- 
tion. Knowledge of the ARP warhead characteristics as well as the target's vul~ 
nerability to each of these effects is essential. Consequently, a preliminary 
analysis is required of the vulnerability of the target to the ARP warhead. 
Fraguentation effects are provided in either of two distinct formats: a P, grid 
which yields conditional kill probability as a function of burst point proximity 
to the target, or a P, vs R (range) function which provides the kill probability 
data as a function of range only; i.e., azimuth characteristics are averaged for 
each range from projectile burst to target. These functions are provided by the 


MAE program. Direct hit and blast effects are estimated from standard target 
vulnerability analysis. 


The overall kill probability for each Monte Carla sample is based upon these 
individual effects and is computed as: , 


Py = 1-C1-Ppg) (L-Pppp) (1-P gis) (1-Pp) 
where Pou = direct hit kill probability, 

Pepg ™ radar blast kill probability, 

PaLsT = vehicle blast kill probability, 


and Py, = fragmentation kill probability. 


Coordinate Syatem 


The simulation uves a rectangular coordinate system whose origin is at 


"ground zero of the target center of vulnerability. Target heading establishes 


the negative range direction (-R); positive deflection (D) is to the left 


(driver's side) of the target; height (H) is measured from the ground (positive 
up) (fig. 2). 
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Figure 2. Coordinate system 


Attack Angle 


The attack angle is the combination of both elevation and azimuth angles 
which define the direction of the incoming ARP with respect to the coordinate 
system for the target. Azimuth is maasured from the negative range direction 
toward the positive deflection. The elevation angle, %, is measured from the 


{ horizontal in the positive height direction (fig. 3)- Azimuth can be either 

| fixed or chosen randomly for each simulation. Elevation is chosen from a 

' Gaussian distribution with a specified mean, -,, and standard deviation, %,. The 
attack angle orients the direction of the ARP flight path (trajectory). 


Figure 3. Attack angle 


Guidance Errors 


Guidance errors are Gaussian and are specified by either the standard devia~ 
tions in deflection and height or CEP in deflection and height. These errors are 
defined in the plane normal to the ARP trajectory and passing through the homing 
point. The location of the guidance plane and the selection of a sample trajec- 
tory through the point (GR, GD, GH) are illustrated in figure 4. The determina~ 
"| tion of the point (GR, GD, GH) is as follows: 
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First, the homing point (GMR, GMD, GMH), defined in the target coordinate 
system (R,, D,, H,), is rotated through the azimuth angle, A. 


GMR' = GMR*cos(A) - GMD*sin(A) 
GMD' = GMD*cos(\) + GMR*sin(A) 


Then GR, GD, and GH are defined based on the sampled errors about the rotated 
homing point. Then 


where H, D ave random normal deviates with 2h = 0, 9= 1, 


GR = GMR' + H*O, *sin(W) 
GD = GND' + D*o, 
and, GH = GMH + H* 0, *cos (“) 


where GR, GD, GH are in the Rp dD, (projectile) coordinate system and Oe o 
are the standard deviations in height and deflection, respectively, of the guii- 
ance error in the guidance plane (Hy » Dy) 


Fuzing 


Six options are available for primary fuzing; both point detonating (PD) and 
proximity (VI) backup fuzes can be considered. Each of the primary fuzes is 
described below: 


Gaussian Fuzing Angle 


Fuze glitter points are specified on the target and a single glitter 
point is selected as either the first glitter point encountered or, optionally, 
chosen randomly for each simulation. When the angle between the flight path and 
a line from the ARP to the glitter point is equal to the fuzing angle, $, the 
point on the trajectory at the vertex of the angle is taken to be the fuzing 
point (fig. 5). The fuze angle for each simulation is selected from a Gaussian 
distribution as, 


od = Hy = NES, 


where V is a random normal deviate with 2b = 0 and Oo = l. 
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Uniform Fuzing Angle 
Identical to the Gaussian fuzing angle except that ¢ is chosen as uni- 
1, formly random between specified limits for each simulation. 
“ Ls Linear Fuzing 
al ; 
pn k 
wk ; Fuzing occurs at some distance along the ARP flight path measured from 
i the guidance plane. The distance along the flight path is chosen from a Gaussian 
. : distribution with a specified mean, ast and standard deviation, o% (fig. 6). 
= Given the ARP terminal velocity, linear fuzing can be used to represent a time 
| fuze where time is measured from the guidance plane. If Hes O, represent the 
an Gaussian parameters for a time fuze, then where Vp is the ARP terminal velocity, 
“ by = Vy*H, and oO = V,*O.. 
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Height Fuzing 


Fuzing occurs at a specific height above the ground. Height is chusen 
from a Gaussian distribution where the mean and standard deviation are 


specified. The point on the ARP flight path which corresponds to the selected 
height is the fuzing point. 


VT Fuze 


A VI fuze functioning distribution is considered by specifying the cumu- 
lative distribution function of fuzing height. A fuzing height is chosen accord- 
ing to sampling from that distribution and the fuzing point is the point on the 
ARP flight path which corresponds to the selected height. 


PD Fuze 


The intersection of the flight path with the ground establishes the PD 
fuzing point. 


All of the above described primary fuze options can have either a PD or 
VI backup fuze. The backup fuze is used if a test for primary fuze functioning 
fails; otherwise, the primary fuze establishes the fuzing point unless a VT back- 


up fuze point occurs at a greater height than the height component of the primary 
fuze point. 


Target 


The physical dimensions of the target are represented by a group (up to 5) of 
rectangular parallelepipeds (fig. 7) with the center of target vulnerability 
located over the origin of the ARP terminal coordinate system (R,, D., H,)- 


Burst Point 


In all cases, once the fuzing point is found, a check is made to ascertain 
whether the target has been penetrated in order to reach that fuze point. If 
such penetration is found, the first penetration point is taken as the warhead 
functioning burst point (in this case, a direct hit burst point). Since the 
burst point is established in the rotated coordinate system (through the azimuth 


component of the attcck angl), prior to determining kill effects, the burst 
point is rotated hack into the target coordinate system.. 
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Figure 7. Target description 
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Direct Hit 


Z£ the burst point of the ARP is found to be at the surface of a parallel- 
epilped representing a face of the target, a direct hit is deemed to have oc- 
curred. 


Blast 
Blast kills can be estimated for both the target vehicle and radar antenna. 
Target Vehicle Blast 


A table of blast radius versus burst height must be provided (fig. 8). 
If the burst point occurs within the radius specified for the determined height 
of burst, then a blast kill of the target vehicle is deemed to have occurred for 
that sample simulation with probability, p (fig. 9 and User Guide, app B)- 


Radar Blast 


A function of the form illustrated in figure 10 must be provided for 
this option. This function defines radar blast kill probability as a function of 
range from the antenna to the burst point. For each simulation, radar blast kill 
ig determined from the specified function. 


Jragmentation 


Fragmentation effects are determined from the results of preliminary MAE 
analysis of the fragmenting warhead. The MAE computer code is described in ref~ 
erences 1 and 2. The MAE program computes conditional kill probabilities as a 
function of burst point proximity to target center, burst height, attack eleva- 
tion angle, and projectile terminal velocity. With the MAE code for a given 
terminal scenario for each of several burst heights, a suitable representation of 
the fragmentation P, function can be described. For each burst height, a P, grid 
is computed which provides the basis for the construction of a P, box grid about 
the target center. It is then a simple matter of interpolating in the range, 
deflection and height directions as well as for elevation angle to estimate the 
fragmentation P, for the actual burst point (fig. 11). Fall-off P, along the 
edges of the Py. box is assumed to be linear out to a specified limit; that is, a 
limit is specified in the range, deflection, and height directions at which the 
fragmentation Py drops to zero. 


It is important to note that the fragmeutation kill probabilities generated 


by the MAE program are based on vulnerability data averaged over all attack azi- 
muths. Also, P's are determined by the MAE code by computation of the proximity 
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Figure 10. Radar blast function 
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of the burst point to the target center of vulnerability poiat. This point-to- 
point relationship is deficient for narrow spray angle munitions in close prox- 
imity tc the target. AJso, since ARPSIM assumes a particular attack azimuth, the 
assumption is made that. for the purposes of the ARPSIM model, the average vul- 


nerability of the target can be used to represent the vulnerability for any par- 
ticular attack azimuth. 


Ag an alternative to the Py grid box, the MAE program can be used to gener- 
ate a P,-versus~range function, where the P, is averaged over all target azi- 
muths. ARPSIM can utilize this function to interpolate for P, based upou the 
range from the burst point to the target center. The Py, ~versus—range function 
can be generated for various burst height and elevation angle combinations. This 
approach is not recommended with directional warheads. 


When the MAE program is used, the blest option available with the MAE code 
should not be used. 


MONTE CARLO ESTIMATES 


The program flow procedures are followed for each simulation to provide 
est::mates for direct hit, body blast, radar biast, and fragmentation effects in 
the form of kill probabilities, P,. Estimates of these kill probabilities are 


computed by using ¥ 
P(t) 
imi * 


PL (n) = n = sample size 


n 


for each of the kill mechanisms. The combined kill probability is computed for 
each sample using’ 


P(A) = LM[1-Pyy(4)] (1-Papg(1)] [Py gy (1) 1 [1-Pp (4) 


These overall kill probabilities are averaged for each individual component 
kill probability to give Monte Carlo estimates of the effectiveness of the indi-~ 


vidual kill mechanisms as well as the overall probability of cefeating the tar- 
get. 


CONCLUSIONS 


The ARPSIM model can be used to provide both weapon designers aud effective- 


neas analysts with an assessment of the potential for the ARP system. As a 
design tool, ARPSIM provides insight into the contributions of guidance and fuz~ 
dag policies to the overall performance of the ARP warhead. ARPSIM does not 
simulate the guidance and control or radiation sensing mechanisms. ARPSIM does 
provide a means to parametrically assess the relative importance of various per- 
formance leveis of the guidance, fuzing, and warhead functions. By providing 
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effectiveneas information for a host of performance capabilities, ARPSIM is a 
useful tool to aid in exploiting those elements of the system which provide the 
greatest payoff in terms of system effectiveness. ARPSIM can also be utilized to 
provide data for systems anaiyses once performance criteria for guidance, fuzing, 
and warhead functioning have been firmiy established by weapon design. 


The following upecific assumptions and limitations are imbedded within the 
ARFSIM model: 


1. Target is engaged in open fiat terrain. 


2. ARP terminal trajectory is linear with the longitudinal axis of the 
projectile collinear with the trajectory. 


3. The target configuration can be adequately represented by an aggre- 
gate of rectangular parallelepipeds. 


4. Fragmentation effects can be estimated with the use of either a P, 
box or a P,.~versus~range function generated by the material lethal area program 
based upon vulnerability data averaged over all azimuths. 


RECOMMENDATIONS 


The computer code follows a sequeiice of steps for each sample simulation. 
Any cf these steps cun be treated as a separate functional module (fig. 1). The 
degree of simulation detail can be changed by developing wore complex modules to 
either increase simulation accu”acy or expand modular function. The consequences 
of either improving the model's resolution or expanding its scope are au increase 
in computer processing time and a resultant increase in the cost of analysis. 
These consequences must be weighed against the advantages to be gained from the 
refinement of the model. 


Some refinements which might be of merit include the direct computation of 
fragmentation effects (rather than use the results of precomputations with the 
MAE code) and the capability to define a complex target array consisting of a 
multiplicity of target elements. 
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This user guide is intended to aid those who have access to the ARRADCOM CDC 
6500/6600 central computing facility via INTERCOM and BATCH wmode processing. 
Others who may wish to use or modify ARPSIM for operation on a different computer 
system should also find this guide informative and helpful. 


For assessment of fragmentation effects with ARPSIM, it is first necessary 
to generate files containing fragmentation P, data as determined by the materiel 
MAE code (ref 1). There are two alternate forms that the MAE-produced PL data 
may take for use by ARPSIM: 


1. AP, grid where grids are defined for the ARP terminal elevation 
attack angle for up to four different burst heights. 


2. A P, versus range table defined for these same terminal conditions. 


For directional fragmentation patterns, the P, grid format provides a better 
estimate of the effects produced by the non~symmetry of the warhead effects pat- 
tern. P, functions produced by the MAE code are developed as follows: 


Py Grid Function 


Several options exist with the MAE code described in reference 1 which allow 
the user to define the bounds of the P, grid in a variety of ways. it is impor- 
tant to note that ARPSIM is limited to a grid size of no more than 20 cells in 
either range or deflection directions. It is quite possible that fragmentation 
effects for an ARP warhead might exist at ranges far in excess of the actual miss 
distance from the target being attacked. For this reason, the user is advised to 
analyze the guidance errors and fuzing scheme being considered in order to deter- 
wine practical limits to the size of the P, grid. Input data for the MAE code 
are often {n units of feet, whereas the P, grid boundaries which are output are 
metric. Also, ARPSIM can be used with any consistent set of units, although it 
is recommended that the metric system be used. It is advisable then to pre- 
determine the practical range for a P, grid and then use an option with the MAE 
code to define the limits of the P, grid. 


When using the OVRLAY code described in reference 1 to make MAE calcula- 
tions, the MTRX option should be called for but not actually used; that is, the 
MIRX input data should consist of a blank card. For the user who is not familiar 
with the OVRLAY system of computer codes, it is a system that was established to 
provide users with the capability to make single computer runs beginning with raw 
fragmentation data continuing through MAE computations and the development of P, 
grids, and culminating with estimates of artillery system effectiveness against 
certain target arrays. The overlay technique is used to combine a number of 
computer codes devoted to these analyses. The MIRX option signals the MAE code 
to produce a Py grid on a file named TAPE4 in forwats waich are compatible with 
both the MTRX and ARPSIM codes. For this reason, the user should call for the 
MTRX option when using the OVRLAY code, and then provide only a blank card as the 
input for the MTRX code. By doing this, the user will normally terminate the 
OVRLAXY code and will have defined a TAPE4 file consisting of a string of P, 
grids, one for each burst height. It is advisable to save the TAPE4 file as a 
permanent file for future recall of the data, as necessary, when using the ARPSIM 
code. 
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If P, grids are being venerated for several (up to three) different attack 
elevations for use by ARPSIM, each elevation angle data set should be generated 
by a separate MAE run. Then, when recalling the P, grid files, define tne data 
on file TAPE2 for the lowest angle data, TAPE3 for the next lowest, and TAPE4 for 
the highest. Burst heights should always be computed in the order of lowest to 
highest. 


For users who do not have access to the MAE code or who will use an 
alternate code to generate P, grids, the files TAPE2, TAPE3, and TAPE4 should 
contain, sequentially, a card image data record in format (213) indicating the 
number of grid coordinates in range and deflection. 


Next, are two card sets in format (10F7.1) where the first set defines the 
range coordinates of the grid boundaries and the second set defines the deflec- 
tion boundaries. Bounderies are defined from lowest to highest values. Follow- 
ing these data sets are the P,‘s associated with the grid in format (10F7.5) 
where P,‘s are given first for the first rauge cell (lowest grid bracket) for 
each of the deflection cells (again, beginning with the lowest bracket) and pro- 
ceeding through all range brackets in the same manner. All P, grids are defined 
this way for each burst height in order of lowest to highest burst height. 


PL Versus Range 


An average P, versus range function (table) can be used if the number of 
‘ranges is no greater than 200. Format for data entry is (F8.3, F8.5) where the 
first item is range (usually in meters) followed by the corresponding average 
P,* The MAE code can generate this table on a file named TAPEL5. These files 
can be saved, like the grid files, and recalled when using the ARPSIM. These 
files when used other than with the MAE code or when recalling MAE-generated 
files, are defined like the grid files, i.e., lowest angle data on TAPE2, next 
lowest on TAPE3, and highest on TAPE4. Each burst height (up to four) has its 
own table defined beginning with the lowest burst height and stored sequentially 
on each file. 


Following definition of the P, functions on TAPE2, TAPE3, and TAPE4 (as 
required), the ARPSIM can be exercised using a teletype (TTY). Preliminary steps 
required to run ARPSIM on the ARRADCOM computer in INTERCOM mode are as follows: 


INTERCOM Mode Setup 


The following sequence is required to access the ARPSIM code and begin 
its execution: 


LOGIN. 

««efollow normal login procedures 
COMMAND ~ ETL,500. 

COMMAND ~ FETCH,ARP,BWEBSTER. 

COMMAND - ATTACH,T,TAPELFILE,ID=“your id. 
COMMAND - COPYBF,T,TAP#1. 

COMMAND - RETURN,T. 
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COMMAND - ATTACH ,TAPE2,... 
COMMAND - ATTACH,TAPE3,... 
COMMAND ~ ATTACH,TAPE4,... 


COMMAND ~ ARP 


The sequence from ATTACH,T... through RETURN,T. is only required if a previously 
defined set of basic inputs is to be used as a basis for this run. Also, the 
sequence ATTACH,TAPE2,... through ATTACH,TAPE4,... is required only in accordance 
with the requirements to estimate fragmentation effects and the diversity of 
attack elevations required. 


In response to the commard ARP, the user will be given the opportunity 
to produce a summary input guide. Following that, the user will be asked whether 
a file named TAPEL is to be used as the basis for input data. This option is 
provided as an aid to the user who expects to make several computer runs with the 
model using the same basic input data set. The ARPSIM code has a built-in input 
editing routine which continually redefines the file TAPE] to be the current 
basic input data set. The user who wishes to make additional runs with a basic 
data set merely has to define the current data set and then, after ARPSIM has 
been run, the TAPE] file is stored on a permanent file for iater use as with the 
ATTACH,T... through RETURN,T. sequence described above. If a basic data set is 
being used, then the initial input conditions are listed. Then, in all cases, 
the user is asked to ENTER DATA OR END - . In response to this command the user 
begins to enter “word” type data to either initialize a data type or change a 
data type. Word type data which can be entered are defined according to general 
function in the section which foiiows. Formats are (A4,F10.4). 


“Word” Type Data 


This section is divided into functional areas as follows: 


Guidance Data 


NGER ,n. NGER signifies the number of puidance error data sets to be 
input. The value of n equals the number of different guidance 
error sets to be analyzed. 


NCEP ,1. If guidance errors are input as standard deviations in both 
deflection and height, omit this set. If errors are input as 
CEP, then include this set. Note that in all cases errors are 
defined in a plane passing through a homing point and normal 
to the ARP flight path. 


Fuzing Scheme 


FZAM,n. FZAM signifies the use of the fuzing angle primary fuze where 
n is the mean value of the fuze half-vertex angle; i.e., n is 
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FZAS ,n.- 


F2ZTM,n. 


FZTS jn. 


PKPF ,n. 


PDVT plie 


GLTR,n. 


OMEG ,n. 


OMGS ,n. 


TGTC,n. 


the mean angle from the ARP trajectory to the fuzing glitter 
point at which fuzing will occur. Units are degrees. 


FZAS signifies the standard deviation of fuzing angle 
associated with the mean value defined by FZAM, where the 
value n is the standard deviation. Units are degrees. 


FZTM signifies the use of # linear (or time) fuze where the 
sign of the value of n indicates whether the fuze operates in 
the vertical direction or along the trajectory. A negative n 
signifies the vertical option. The value of n is the mean 
distance from the guidance plane (or initial fuzing point if 
used in conjurction with the FZAM option) in the negative 
range direction where fuzing occurs. With the vertical 
option, the distance is measured from the ground. A time fuze 
operating along the ARP trajectory can be simulated b. con- 
verting the values to distances by using the known AR: termi- 
nal velocity. 


FZTS defines the standard deviation associated with the FZTM 
data in all modes. 


The value of n is the probability that the primary fuze 
(options described above) will function. 


Selects the backup fuze option. The value for n is O for a PD 
(ground burst) backup and is the number of entries in a height 
versus probability table (up to 5 values) to define the VT 
fuze functioning distribution. 


Specifies the glitter points uged by the angular fuzing 
function option. If nun is 0, the fuze functions relative to 
the point (0,0,TGTC) where TGTC is the center of target 
vulnerability. If n is non-zero, the fuze functions relative 
to one of the n input glitter points. A positive u signifies 
that the fuzing glitter point is selected randomly; a negative 
n signifies that the first glitter point encountered will 
cause fuzing. 


Terminal Conditions 


The elevation angle measured from the ground is chosen from a 
normal distribution with mean value n. 


The standard deviation associated with OMEG is input as n. 


The center of target vulnerability is input as a height above 
the origin at (0,0,n). If direct hit effects are not being 
analyzed (direct hit boxes are not defined), then the vehicle 
blast effects are determined based on the range from the burst 
point to (0,0,TGTC). 
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Eo 8 DHAZ,n. The azimuth angle-of-attack is n and is measured from the 
oo negative range axis in the direction of the positive 


ie deflection axis. Units are degrees. To choose the azimuth 
ee ke uniformly random between O and 360 degrees, set n = -l. 
4 DUDR,n. The dud rate of ARP projectiles is given as n, where a 5% dud 


rate corresponds to n = 0.05. 


General Conditions 


SAMP,n. The number of Monte Carlo samples is n. 


oot PRNT,1. Specifies that only a final summary of results is to be 
output. 


SRNG,n. Tables of average combined P, can be output as a function of 
azimuth, elevation and range as well as averaged over non~zero 
results obtained in the angular bins for each range. The 
value for n is the upper limit (defaults to 100) for range 

L information. The range scale is logarithmic and includes 10 

ed bins, beginning with the minimum range obtainable (considering 

(ae direct hit implications) and ending at n. 


Fragmentation Effects 


PKNH,n. Specifies the number of heights, n, at which fragmentation 
effects are provided (either as PL grids or P. versus range 
tables). Must not exceed 4. 
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PENA,n. Specifies the number of elevations, n, for which fragmentation 
effects are provided. Must not exceed 3. 
For n = 1, effects are on TAPE2. 
For n = 2, effects are on TAPE2 for lowest angle data and on 
TAPE3 for highest angle data. 
For n = 3, effects are on TAPE2 for lowest angle data, TAPE3 
for middle angle data, and TAPE4 for highest angle data. 


xf See ere 


sTageenaere 


‘ 
| 
4 
4 


7S 


FUNC,1. Selects option to use P, versus range tables for fragmentation 


ef effects in place of the P, grids. 
io | Direct Hit Effects 
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DHIT,n. specifies the number of target boxes to be input tc 
approximate the shape of the target for purposes of computing 
direct hit effects. Boxes are defined relative to (0,0,0) and 
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PKDH,n. Direct hit Py. if a direct hit is achieved. If n~ 0, P, is 
defaulted to one. 


Blast Effects 


PKBL,n. Specifies the blast P, if the burst point is within a range 
specified by the BLST data of the surface of any direct hit 
box. If direct hit boxes are not used, then range is 
calculated to the point (0,0,TGTC). 


BLST,n. Specifies the range from the divect hit surfaces or the point 
(0,0,TGTC) within which the blast PL against the vehicle body 
is that given by the PKBL data. ‘To ener a table of blast 
ranges versus burst height, enter a negative value for n which 
corresponds to the number of entries in the blast range versus 
height table (may not exceed 5). 


RADR,1. Include to compute blast effects against radar antenna 
separately from vehicle blast. 


End of Word Data 


END Must always be included at the end of the “word"=type data 
entries. 


After all “word"-type data have been entered, the code will ask for 
certain data which are required by some of the options chosen by the “word” 
cards. These additional input requirements are discussed in the following sec— 
tion. All data are free~formatted. 


Guidance Data 


Either pairs of deflection and height standard deviations are entered 
or, if NCEP,1. data is entered in the “word” section, then the guidance errors 
are input as CEP's. 


The homing point coordinates follow the guidance error inputs. The 
homing poiut is generally the coordinates of the center of the radar antenna. 


YLlrect Hit Boxes 


The limits of the dimensions of each direct hit box are input for range, 
deflection, and height, respectively. For example, for a direct hit box centered 
at the origin and having a length of 20 meters, a width of 10 meters, and a 
height of 5 meters, this data would be input as ~10,10,-5,5,0,5. 


Radar Data 


Radar antenna coordinates are entered for the purposes of radar blast P, 
computacion. 


Following the entry of the radar coordinates, values are entered for two 
ranges, Rl and R2, which define the radar blast P,. function as being one out to 
Rl and declining linearly to zero at R2. 


Fragmentation 


Heights are entered beginning with the lowest value and corresponding to 
the burst heights used for the MAE computations. An additional height is input 
last and corresponds to that height at which all fragmentation P's are zero. 


Following the height datu, two values are input corresponding to the 
distances beyond the edge of the P, grids where the fragmentation P, becomes zero 
in range and deflection, respectively. 


Elevation angles are entered next, beginning with the lowest angle and 
corresponding to the angles for which the MAE code was run tv produce the frag~ 
mentation P, data. 


VI Backup Fuzing 


A table of probability of fuze functioning at helght less than or equal 
to height, H, is used to generate VI fuzing data. Up to five heights are input 
followed by probabilities corresponding to the: probability of fuze functioning 


between the respective height and the next lower height. Ideally, probability 
values should sum to unity. 


Glitter Points 


Glitter point coordinates are entered for each glitter point. All coor~- 
dinates are relative to (0,0,0) of the target. 


Blast Data (Vehicle) 


If the blast~distance-versus-burst-height cption is chosen (negative nu 
on BLST,n data), then n pairs of blast distance, height are entered. 


This concludes the input requirements for using the ARPSIM model. Word type 
data can be changed or input in any order. Required additional data will be 
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prompted from the user by the code. The user is always given the option of list-: 
exception of the fragmentation P. data) or 


ing the current data set (with the 


changing the data set prior to actual computations. When the computations are 


completed for all cases, the user is given the opportunity to run additional 
cases based on the current data sets. 
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The following example, provided as a supplement to the User Guide 
type of material generated for a typical ARPSIM run: 


dix A, denotes the 


SESARATTASTALSLALATPSSATSSSSESALESTSSTSSSSTIGI SY 
: ANTI-RADIATION SIMULATION PROGRAM ~ 971/780" 


RESEATATELESTTAESATRAAASSSTSSASRTSEALSSESSIAALSS 


ee NOTES ALL COORDINATES ANE DEFINED RELATIVE TO" 
: ORIGIN AT GROUND ZERO OF TARGET. ° 

$ COORDINATE SYSTEM 15 RECTANGULAR, ° 

© TARGET HEADING 1S NEGATIVE RANGE, ° 

ty DRIVER SIDE (L) 18 POSITIVE DEFLECTION, * 
. HEIGaT 15 MEASURED FROM GROUND. °® 


RSEAUTKSTASLASSSATAAKAAKALTAAG SALTS L TEATS AIST 


DATE - 08427780 
TIME ~ 33.47.13. 


KESASTUAESRRSCRISSATCIULTTARAALT SASS LSRLSIITS AT 
*DO YOU WANT A LISTING OF CODE NAMES? °¥ 


SSTALSTSLTSESEATTISELSTAT TAL TASSTALBATLCLESTTLSAT 


“OMEG - MEAN ATTACK ANGLE? 
*ONGS - ATTACK ANGLE STO DEV* 
* NOTE’ FOLLOUING GUIDANCE ERROR PARANETERS® 
° ($IGD,SIGH) ARE MEASURED®* 
. IN PLANE NORMAL TO TRAJECTORY AND® 
. PASSING THROUGH HOMING POINT® 
NGER - NUMBER OF GUIDANCE ERRORS TO CONSIDER* 
ENTER HOMING POINT (R,D,H), GUIDANCE® 
° ERRORS ARE DISTRIBUTED ABOUT HOPING PT.°® 
*NCEP - 1., IF CEP IS INPUT FOR GUIDANCE ERROR SIGNAS® 
SFUNC - de, IF OPTION TO USE PK VS. RANE DA“AGE® 
IN PLACE OF PK 8OX FUNCTION IS SELECTED‘ 
YOU MUST DEFINE PK US R DATA FOR EACH’ 
HEIGHT LAYER SPECIFIED BY PKNK CARD’ 
AND EACH ANGLE SPECIFIED bY* 
PKNA CARD. ® 
FZAN,FZAS,FZTM,F2TS - FUZING ERROR OPTIONS® 
FZAM - MEAN ANGLE AT UKICH FUZING OCCURS ON® 
INTERCEPT’ 
F2AS - STD DEV ASSOCIATED WITH FZAM* 
NOTE: FOR UNIFORM FUZING ANGLE BETWEEN FZAM® 
AND FZAS, ENTER R& NEGATIVE VALVE FOR FOAR® 
FUZE ANGLE WILL BE CHOSEN UNIFORALY RANDOM® 
BETWEEN POSITIVE FZAM AND FZAS* 
NOTE! FUZING PLANE PASSES THROUGH FLZING GLITTER® 
POINT NORMAL TO SAMPLE TRAJECTORY® ; 
F2ZTM - MEAN DISTANCE FROM GUIDANCE PLANE AT UHICH® 
FUZING WILL OCCUR ALONG TRAJECTORY’ 
*NOTES ENTER A NEGATIVE FZTN FOR HEIGHT FUZING ° 
UITH MEAN HEIGHT ABS(FZTN) * 
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* F2T8 : oY arb ade ae WITH FZTR* 
*SRAP - BP LE 

RUMDER OF HEIGHTS AT UHICH FRAGNENTATION® 
PK DATA UILL BE DEFINED* 
3 MOTE! PkKNn ¢ S° 
*PKNA ~ NUNBEW GF ELEVATION ANGLES FRAGMENTATION® 
° PK DATA UILL BE DEFINED FOR® 
$ NOTE! PKHA ¢ 4° 
*PKPF ~ PROBABILITY OF PRIMARY FUZE FUNCTIONING* 
“PDUT = ©, FOR PD BACKUP, NYT FOR UT BACKUP FUZE® 
: UHERE NUT © NUMBER OF UT BURST HEIGHTS’ : 
*GLIR ~ ©. IF PRIMARY FUZE FUNCTIONS RELATIVE TO® 
. CENTER OF TARGET, NGLT IF PRIMARY FUZE® 


¢ FUNCTIONS RELATIVE TO ARY ONE OF NELT® 
¢ EQUALLY LIKELY GLITTER POINTS .° 

. SET NGLT NEGATIVE TO PICK FIRST: 

. POINT ENCOUNTERED. * 


*6RNG = MAXIMUN RANGE FOR COMPUTING PK VS RANGE* 
Heed 1. TQ PRINT SUMMARY ONLY, @. OTHERVISE* 
&. TO PRINTOUT PROGRAM DESLGGING DATA’ 


ee 
ay 
ax 
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: DBUG © 1, GUIDANCE & FUZING £ \Ta® 

. DbuG © 2, DIRECT HIT PENETRATION DaTas 
. 08UG © 3, HONING ANGLE DaTa‘ 

b OBUG = 4, PK BOX DATA* 

: O8uG + &, PK GRIDS’ 

: BBUG © 6G, PK US R DATA® 


*TGTC = HKIGHT OF TARGET CENTER ABOVE GROUND® 
*DUDR = DUD RATE OF PROJECTILE, EXPRESSED 45 & FRACTION? 
*OAST + DIRECT MIT OPTION, HUMBER OF TARGET BOKES® 


. IF DHIT 2S OMITTED waND BLSY 15 KNCLUDED, * 
7 BLSY IS RADIUS ©2I0N (@,0,TCYC) UYMIN: 
° GHICH PKBLST 2 4.° 


SPKDH = DIRECT MIT PK (MH, © 4,2° 
“PKBUL =~ ULAST PK (Q@, © 4,29 


*RADR = 1., DEFINE FUNC FOR BLAST KILL OF RADAR ONLY® 
U AND READ IN RADAR ANTENNA COORDINATES, * 

: TO DEFINE FUNC, SPECIFY RI AND RO, 

. WHERE BLAST PK I$ 1 GUT TO R31 anD’ 

- CECLINES LINEARLY TO @ at Re.° 


*DHAZ = AZIMUTH ANGLE OF ATTACK OFF FRONT OF TARGET® 
: TuJHRO ORIVER SIDE. SET TO -3. FOR RAtIDOM:® 
*BLST = ULAST RADIUS WITHIN UHICH VEHICLE PePKBL® 
"MOYES TD ENTER BLAST RADIY US, BURST HETAHT, ° 

: ENTER NEGATIVE NUMBER OF BLAST ,HGT UATRS® 


r TH PLACE OF VALUE UF BLST. PINS OFF 
: BLAST, WGT ARE ENTERED IN ASCENDING ORDER® 
i OF HEICHT.® 


“COORDINATE SYSTEM 15 RECTANGULAR, ° 

“TARGET HEADING 1S NEGATIVE RANGE, ° 

“DRIVER SIDE (LEFT) IS POSITIVE DEFLECTION, ° 
*MELGHT IS MEASURED FROM GROUNK® 

“ENTER DATA BY ENTERING COLE NaAnEet 

“FOLLOWED BY @ COMMA AND THE VALUE IN FLOATING® 
“POINT FORMAT. YO END DATA ENTRY, ENTER ° 

“THE WORD END IN COLUMNS 1-3° 

Ch YQU WISH TO INITIALIZE DATA FROM SAVED’ 
“DFTA FILE (TAPEL)2 *¥ 


*IMITIAL INPUTS © * 
FZAaN 


70.0006 
PROUOK 1.860 
PKBL 1.998 
FZAS 1¢.000 
OMG 10.008 
NGOER 3,080 
NCEP 1.300 
FUNG 1,080 
DHIT 2.920 
SANP 160.080 
PKH 4.008 
PKNA 2.880 
PDVT 5.820 
PKPF -950 
GLTR 3.028 
SRNG 10@.02e 
TGOTC 16.000 
UUDR -05¢ 
BLST 3.008 
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°00 VOU WANT 76 CHANGE ANY DATA? = +y 
‘ENTER DATA OR END - “RBDR, 1. 

‘ENTER DATA OR END - “END 

ENFER RAGAG “ANTENNA <CSONDINAYES QR, B>Wa- @eUATIUL 
YU TARGET GROUND ZENO. -@.,0..20, ° 
ENTEW R1,R2, UMERE RADAR BLAST PKey 
OUT TOW} AND LECLINES LINEARLY 

TO ZERO AT RZ ~18.,20, 

"00 YOU UAKT CURRENT DATA LISTED? °N 


"BO YOU VANT TO CHANGE ANY LATA? = KN 


+OdSe NSLAP 


BETA" 


ee rer) 


60.0 
14.0 
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SSACOARE ST CAAUGKEENELSSESISSAVSSASAUSVSSS SITS Hayes 
FINO We TILTS* 


RUVROGLSSASHERKERARETLERSIIALAIKARN SEXTET TEESE 
°OO YOU WANT PK US R, a@LPHA, BETA? ty 


75.0 
39.0 
45.0 
60.8 
o8.0 
25.6 
75.0 
4$.0 
60.0 
45.8 
60.0 
49.0 
75.0 
GQ.6 
99.8 
68.8 
69.9 
75.8 
62.0 
8G.0 
45.¢ 
u0.0 
275.0 
45.0 
62.0 
see 
eed 
75.0 
45.0 
bo.0 
45.0 
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Toe 
Ae ot 
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‘3% = a BEES SN, ET TREE SE EE a ke 
ay =~ - 
: alee a ee bt aeee 
aes se ee eam =f 
. J i 
an | 
‘ 
KY i 
Be | 
ee SPPESTISGresiressreritttiitiittittitirititt Tritt tt 
: "AUG BK US, Re 
te i Bene nenacuwaat . 
n: 1.0008 11.8 
: +9026 12,7 
ye. , Pas Tn 14.4 
; VIse 16.2 
- .33e 1$.7 
ay =A 62454 25,1 
: eh hds 32.5 
O15 40.7 
SICTSITSLSTCHLACSCUSSERTSSVISRELTSTSIEEKINVTVHSISIALE 
iy. 
. “PINAL RESULTS* 
PK «© ,@978 PKSD = .6409 NSAMP « $60 
CRESTSRSSALTRASTSTSLESASANSTACEAECRNS CREATES T ETE 
“PO YOU WANT PK US R, ALPHA, BETH? *¥v 
+ PK R ALPHA beta: 
\ eae wtws wan te Sunreeruseenwraseasn @eu2eeeeue =ewertaee ° 
1.002@ 31.8 120.0 - 150.0 16.0 - 30.9 
: 1.0020 31.8 120,09 - 150.0 30.0 + 45.a 
1.8028 11H 120.8 - 158.0 45.0 = 60.0 
1.0000 11.8 220.0 - {5¢.0 6e.w = 95.9 
1.0000 11.8 15,0 - 1hd.e $2.8 6 45.9 
‘ 1.0086 11.8 15@.@- 1hd.¥ 45.0 = Gu.d 
‘ 1.0000 11.8 15@.@- if0.0 Be. = 75.8 
1.0000 41.8 218.0 - 240.0 ON cr De) 
1.0000 12.7 {$2@.0 = 150.0 Bd = aha 
ail 16573) 12.7) 1SR.0 = 1RA.8 15.0 5 JH.8 
1.0000 12.7 15e.0 - 180.0 69.0 - 24.0 
iy 1,8600 12.7 aid.d- advc.e 45.8 = 66,0 : 
a | 1.0200 9 14.1 Ve = 120.8 EQ. - 75.8 
1.0008 14.4 1288.0 © 150.0 45.8 = 66,0 
he 1.8000 i4.f 128.0 - 159.0 Tr nr are 
,8977 14.4 182.0 = thaee 45.8 = Ge.d 
1.2000 oo 14.$ 150.8 = 18a.e 68.9 = 75.9 : 
L.8Cud 14d 21d,0 = Bence 45.8 = 60.2 
L.se0d «16.2 92.0 2 LEOLY be.e@ = 74.9 
19949) 16.62) eM -  1SOLe 15.2 = bead 
26972 1G.2 152.0 = 1bO.e Wee = aoav 
19549) 16.2 158.0 - 1kd.e | AB e 6OL0 | 
1.0000 16.2 158.0 6 Lae.e 606d = TE, i 
64425-1662 RB = ata.a S.e - 30.9 ny 
$8000 16.2 2id.e - 240.0 45.2 =  6¥.9 Dey 
Lacy = $6.20 eid.e - 84002 60.@- 75.0 1 | 
£9039 15.7) «188.0 - 150.8 45.0 = 6e.0 
1.98000 15.7 1288.0 - 180.0 60. = Ih.a | 
1.0000 15.7 $28.0 + 180.a "5.0 2 99.0 
12305 19.7? 158.0 2 ThMLE iSud © mee 
69559) 15.7 Sea Wee abae 
+BSEG 0 1S.7 = $58.0 - 1B0.6 45.0 - 62,0 ; 
1.8080 915.7 $5d.8 - Fea.e 6I.E - 75,0 
L.9CKO 1567 Ge eae Ge = YEG 
te Ce Ge thd 
»I5ES 1>.? Q1a.d -  Baa.e 45.0 - 6a.8 
POON 15.7 Bre. = Z4aee OI = 75.0 
1.806 2.5 $20.8 - 18g ba.e@ = 75.3 
Oded 25.50 1B thea “GS. = OMG 
vedo OS.) 1G2.8 - eee 19.0 = 30.0 
9414 25.$ 0 tSa.u Sede 45.60 2 Eid 
12338 25.5 LEO de Et eee Qe = 15,0 
1777) BELL SRO Bt OLY 15.0 = 30.0 
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1.000@ 11.8 820.0 = 150.0 15.0 - 38.0 

: 1.0008 11.8 128.0 - 150.0 30.8 - 45,0 

i. 1.0009 11.8 120.0 = 150.0 45.0 - 60.0 

; 1.0008 11.8 150.0 - 180.2 45.2 = 60.0 

Py 16204 «14.8 210.0 - 240.0 16.98 - 30.@ 

: 1.000@ 11.8 210.0 = 240.0 30.8 - 45.0 

! 1.38009 12.7 129.0 = 150.2 30.4 « 45.¢@ 

: 1.0000 12.7 120.@ - 150.0 45.0 - 60.0 

+$632 12.7 158.0 - 198.0 @.2 > 15.8 

©7695 12.7 150.0 - 180.9 30.8 = 45.0 

8241 12.7 150,80 - 182.0 45.9 = 628.0 

P 1.0000 12.7 218.@ - 240.0 45.0 - 60.0 

64675 t4.i 30.8 =. 69.2 17.0 =: 33.@ 

x a 1.0000 14.5 120.0 = 158.8 45.9 - 62.0 

i 66169 14.1 150.0 - 1828.0 15.0 = 328.4 

1.0088 4.1 210.@ - 240.0 45.0 = 63.0 

-98)8 16.2 120.0 = 150.0 45.9 - 66.0 

1.0000 18.2 120.0 - 152.0 60.0 - 75.@ 

14140 16.2 «6158.82 ~ 180.0 15.0 > 30.6 

+6730 «16.2 160.8 - 188.0 3Q.0 = 45.0 

.8445 16.2 150.8 - 139.0 45.0 - 60.8 

1.0020 16.2 150.0 - 1823.0 60.0 = 75.8 

1.00800 16.2 210.8 - 2420.0 45.9 - 62,2 

4.0083 16.2 218.2 - 2428.0 60.0 - 75.0 

19439 19.7 128,08 - 150.0 45.0 = 60.8 

1.0000 19.7 120.0 - 153.0 60.0 - 75.0 

£.0082 19.7 120.0 - 158.0 75.8 > 90.8 

64614 19.7 158.02 = 188.0 15.@ - 30.8 

239359 «19.7 188,02 = 188.0 30.0 - 45.8 

-9249 19.7 150.0 - 182.0 45.0 = GA.e 

1.0002 19.7 158.8 - 180.0 60.0 - 75.¢ 

1.0000 19.7 150.8 - 180.0 75.8 - 90.0 

; 1.0080 °19.7 210.0 - 240,28 60.0 = 75.2e 

i 1.0000 25.1 120.@ = 158.0 6@.0 - 75.8 

: £.0088 5.1 $2e0.@ - 158.0 75.0 - 90.¢ 

i 22502 = 25.1 190.8 - 180.8 16.0 - 30.¢ 

-$088 25.1 1560.0 - 180.9 60.8 - 75.¢ 

: e247 25.1 180,88 - 210.0 0.8 © 15.8 

1878 0 ©62S.1 188.2 - 210.0 15.0 - 30.0 

i ©9106 25.1 188.0 - 210.8 60.0 - 5,0 

1.0082 25.1 210.0 = 240.9 66.2 - 75.a 

29755 33.5 120.8 - 158.28 F1.0 - 75.0 

1173 33.S 1580.0 - 180.8 B.0- 15.0 

“£3916 33.5 158.8 - 180.0 15.0 - 39.¢ 

.2503 33.5 159.0 - 186.2 45.@ - G06 

65454 33.5 150.0 - Ihde 60.8 - 5.0 

$345 33.5 188.0 - 218.28 60.2 = 75.8 

4.0009 33.5 219.8 - 24¥.0 0.2 = 75.0 

| : 194990 46,7 159.8 - 1be.e d.e - 15.0 

68450-4667 189.8 - 218.8 @.c - 15.0 

18389 7.5 150.0 = 180.2 8.82 - 35.9 

i 10029 0-67.50 188.8 - 210.2 @.0 - 15.0 

.8909 180.8 153.0 - 189.28 ee - 15.8 
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RESULTS FOR FOLLOUING CONDITIONS - 


tTen MEAN STD DEV 
ELEVATION 16.0000 6.900@ 
FUZE ANGLE 70,0008 18,0000 
LINEAR FUZE 0. Unde O.du0¥ 
ACT VU TH BG 


SAMPLE SIZE - ine 
HOMING POINY COORDINATES (R,D,H) © @.8, @.8, 10. 


ERROR DATA eK PKFRAG PKRADR PKDHIT PKBLST 
CEP - 3.8 ©7816) 64173) 6.6809) 6708) 6S S00 
CEP = 6.4 6579) 65299) 6 S517. ASO 64000 
SEP = 9.@ 1698) .2613 .5725) .elou w ded 
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A description of the material produced by this particular ARPSIM run 
follows: 


Header information is printed, including the time and date of the run. The 
user is asked whether a listing of input code names is desired (as an aid to 
generating a proper set of inputs). In this example, the code names are printed. 
Next, the user is given the option of starting with a previously developed set of 
inputs which can be changed by a built-in input editing routine. That option is 
invoked for this example. Note that a file named TAPE] must be defined which 
contains this data prior to running ARPSIM. A listing of initial data conditions 
is provided next. The user is then asked whether any data changes are required. 


In this example the user desires to add the capability to estimate radar 
blast effects. Note that only changed data need be entered at this point. The 
code then asks for additional information required by the added data. Having 
fulfilled the data requirements, the user is given the option of listing the 
entire data set again. Following this, the user is given the option of making 
any additional changes or corrections to the data set. In this example no addi- 
tional changes are requested. 


Before proceeding with the discussion of the ARPSIM results for this case, a 
brief run-through is given of the input data set. The FZAM data specifies a fuze 
angle option with a mean value of 70 degrees for the fuze angle. The FZAS code 
specifies a l0-degree standard deviation for the fuze angle from simulation to 
simulation. The PKDH and PKBL data indicate direct hit and vehicle blast P.'8, 
respectively. Attack elevation of 10 degrees is specified by the OMEG card. 
NGER indicates three different sets of guidance errors will be analyzed, and NCEP 
indicates that guidance errors will be input as CEP. FUNC specifies that the 
fragmentation P's will be estimated from interpolations in a set of P, versus 
range tables generated by the MAE code for a combination of burst helght and 
elevation angles. 


Up tc three elevation angle sets can be provided on files TAPE2, TAPE3, and 
TAPE4. If only a single elevation angle data set is provided, then only TAPE2 is 
required. Two elevation angles require both TAPE2 and TAPE3. Each file contains 
P, versus range for identical burst heights, beginning with the lowest burst 
height. That is, if four burst heights have been analyzed by the MAE code (the 
maximum allowable by ARPSIM), each file will contain four P. versus range tables, 
one for each burst height beginning with the lowest height and progressing to the 
highest. i 

! 


In this example, four burst heights were considered for each of three angles 


of fall (elevation angles) as specified by the PKNH and PKNA codes, respectively. — 


SAMP provides the number of simulations to run for each case. PDVT specifies 
that a VI backup fuze is being considered where the height of burst distribution 
for the backup fuze wiil be typified at five burst heights. PKPF specifies that 
the probability that the primary fuze functions is 0.95. GLTR specifies that 
three glitter points for primary fuzing exist. SRNG gives the maximum range for 
a Py, versus range table to be generated based upon the results of the ARPSIM 
run. TGTC provides that the center of target vulnerability is located at 10 (in 
this case meters) above the target origin (0,0,0). DUDR specifies a projectile 
dud rate of 5%. BLST provides a blast radius from the TGTC point within which 
the P, for vehicle blast effect is as stated on the PKBL data above. 


one FRE Es il i dal 


The END code signifies the end of the word type data. The numbers 3., 6:, 
and 9. specify the guidance error CEP's. Following this are the homing point 
coordinates (0,0,10), and the limits in range, deflection and height of the two 
direct hit target description boxes. Burst heights and angles of fall (eleva- 
tions) utilized by the MAE code in generating the P, versus range tables are 
specified next. Then the heights snd probabilities associated with the backup 
fuzing function are listed. Finally, glitter point coordinates are specified. 


Final results are given as the combined kill probability, the standard devi- 
ation of kili probability and the sample size upon which these numbers are based. 
The user is given the option of listing the generated hemispheric distribution of 
computed combined P,'s, where the angle alpha denotes azimuth and beta denotes 
elevation from the burst point to (0,0,0). The range specified is also the range 
from the burst point to the origin (0,0,0). These hemispheric data (only the 
positive elevation angles are considered since negative angles would imply a 
burst belew ground) are averaged over all angular bins for which burst points 
were analyzed to provide a table of average Py versus range. 


The final results are repeated for each case and followed by a summary of 
the results for each type P, considered together with the corresponding error 
data for that case. 


After all results have been given for all cases specified, the user is given 
the opportunity to run additional cases, based upon the same data set. In all 
cases, the contents of the file TAPE1 are always the last data set considered. 
Consequently, if the user wishes to make additional runs with ARPSIM at a later 
time using the same basic data set, then after the current runs with ARPSIM are 
finished, the file TAPE1 can be saved as a starting point for future runs. 


TAPE] can be retained as a permanent file. However, for access at a later 
date, this TAPE] must be attached with a different local file name. Then this 
local file name is copied to a new file named TAPE]. These steps are necessary 
because the ARPSIM code changes the contents of the file [APEl. 


APPENDIX C 


FORTRAN LISTING 


The following FORTRAN listing is 
subject to changes as dictated by 
improvements or modifications to 
the ARPSIM model. 
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